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Abstract. The particles of a classical relativistic gas are supposed to move under 
the influence of a quasilinear (in the particle four- momenta) , self-interacting force 
inbetween elastic, binary collisions. This force which is completely fixed by the 
equilibrium conditions of the gas, gives rise to an effective viscous pressure on the 
fluid phenomenological level. Earlier results concerning the possibility of accelerated 
expansion of the universe due to cosmological particle production are reinterpreted. A 
phenomenon such as power law inflation may be traced back to specific self-interacting 
forces keeping the particles of a gas universe in states of generalized equilibrium. 
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1. Introduction 

In its most elementary version the kinetic theory of a simple relativistic gas relies 
on the concept of N pointlike particles which may interact through elastic, binary 
collisions. Inbetween the collisions which are assumed to establish an (approximate) 
local or global equilibrium of the system the particles move on geodesies of either a given 
spacetime or a spacetime selfconsistently determined by the gas particles themselves. 
Technically, the gas particles are described by an invariant one-particle distribution 
function governed by Boltzmann's equation. The macroscopic fluid dynamics for such 
a system may be obtained in terms of the first and second moments of the (close-to- 
equilibrium) distribution function. An interesting area of application is the dynamics of 
a cosmological fluid. While hydrodynamic concepts such as energy density and pressure 
are assumed to be useful back to times of the order of the Planck time it is less clear 
whether the hot and dense early universe is accessible to a kinetic description as well. 
A gas, however, is the only system for which the correspondence between microscopic 
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variables, governed by a distribution function, and phenomenological fluid quantities is 
sufficiently well understood. This property makes gas universes interesting toy models. 
There exists a considerable body of literature refering to model universes for which a 
kinetic approach is assumed to be applicable (see, e.g.,|Q 0, ^ and references therein). 
One hopes that such kind of approach is also able to give an idea of the relevant physics 
in our real Universe. 

The behaviour of structureless particles moving on geodesies inbetween elastic, 
binary collisions, the latter being subject to the hypothesis of "molecular chaos", 
equivalent to neglegting interparticle correlations, is by now well understood. 
Sophisticated solution techniques have been developed and applied to numerous 
physically relevant situations 0, ||, H, 0, H, H. The usual procedure here is first to 
characterize equilibrium states, i.e., states with vanishing entropy production and to 
relate the parameters of the corresponding distribution function to macroscopic (perfect) 
fiuid quantities. Nonequilibrium situations are then, in a second step, taken into account 
as deviations from equilibrium. 

The concept of a simple relativistic gas has been refined in various ways. For 



example, the gas particles may be supposed to have internal structure [jTO|, |TT|. Taking 
into account interparticle correlations in the simplest possible manner is equivalent 
to an additional stochastic term in Boltzmann's equation, leading to fiuctuating 
hydrodynamics [|l^, |13|. In order to investigate cosmological particle production 
processes, supposedly of quantum origin, "source" terms have been introduced into 
Boltzmann's equation giving rise to corresponding sources in the macroscopic balance 
equations for the extensive thermodynamic quantities (particle number, energy, entropy) 
10, ^ [I^. Under certain conditions these macroscopic sources are equivalent to effective 
viscous fluid pressures. 

Any of these modiflcations and extensions of the kinetic theory of a simple gas is 
supposed to enlarge its ranges of applicabihty and credibihty since it drops at least 
one of its simplifying assumptions. 

In the present paper we investigate a system of point particles interacting through 
elastic, binary collisions but not necessarily moving on geodesies inbetween the collisions. 
Geodesic motion of particles is a highly idealized case. In reality, particle worldlines are 
supposed to deviate from geodesies since the particles will be subject to additional 
interactions in general. We assume here that these interactions may be modelled as 
effective forces on the particles. The kinetic theory for particles under the influence of 
various forces was considered, e.g., in p, ||, |15|, |16|. We will focus on equilibrium states 
of the gas in a simple force fleld. It is well known that a Lorentz force which influences 
the particle motion inbetween the collisions will change, in general, the conditions under 
which the system of particles is in equilibrium, but not the functional structure of the 
equilibrium distribution function. Here we are interested in forces which are quasilinear 
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(in a sense described below) in the particle four- moment a. We will find the most general 
quasilinear force compatible with the equilibrium distribution function of a classical 
Maxwell-Boltzmann gas. The spacetime structure of this force is not assumed to be 
given but will be determined from the (generalized) equilibrium conditions for the 
gas. One part of the force will turn out to be proportional to the fluid expansion 
of the system, another one proportional to the spatial gradient of fi/T where n is the 
chemical potential and T is the fluid temperature. It is characteristic for our approach 
that the force depends both on the microscopic particle momenta and on macroscopic 
fluid quantities which characterize the A^-particle system as a whole. Consequently, 
it describes a selfinteraction of the gas. Physically, the concept of a selfinteracting 
force may be regarded as a mapping of parts of the interactions in a many-particle 
system onto an effective one-particle quantity. We expect this concept to be useful 
in circumventing some of the general problems inherent in attempts to formulate a 
relativistic statistics for interacting many-particle systems (see, e.g., [|l7l and references 
therein). The selfinteracting force neither preserves the particle number nor the energy 
momentum and it will give rise to entropy production. The corresponding source term 
in the energy momentum balance may consistently be mapped on the effective viscous 
pressure of a conserved, imperfect fluid energy momentum tensor. This deepens our 
understanding of the phenomenological description of microscopic processes in terms 
of effective bulk pressures. We will establish the correspondence to previous results 
0, 1^, Q on kinetic theory and particle production which are recovered here in a 
physically more transparent way and with less assumptions. We show that the effective 
viscous pressure which may considerably modify the dynamics of a Friedmann-Lemaitre- 
Robert son- Walker (FLRW) universe, including the possibility of power law inflation, is 
related to specific selfinteracting forces on the microscopic gas particles. Moreover, for 
interesting limiting cases the newly introduced force concept allows us to integrate the 
microscopic equations of motion explicitly and thus to establish an exactly solvable 
model of a selfinteracting gas universe. 

The plan of the paper is as follows. In section 2 we establish the general kinetic 
theory of a relativistic gas under the influence of an arbitrary four-force and recall 
the concept of "generalized equilibrium" . The corresponding equilibrium conditions for 
quasilinear forces are obtained in section 3. Section 4 is devoted to the fluid dynamics of 
a selfinteracting gas in generalized equilibrium. Earlier results on the role of cosmological 
particle production in a phenomenological context are reinterpreted and effective viscous 
pressures are related to forces on the microscopic gas particles. Section 5 gives a brief 
summary of the paper. Units have been chosen so that c = kB = h = 1. 
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2. General kinetic theory 

The one-particle distribution function / = f {x,p) for relativistic gas particles under 
the influence of a four-force = (x, p) obeys the Boltzmann equation [§|, §|, 

PV«-rLpV^ + mF^^ = C[/] , (1) 

where / (x, p) p^rikdY^dP is the number of particles whose world lines intersect the 
hypersurface element n^dE around x, having four-momenta in the range (p, p + dp); 
z, I ... = 0, 1, 2, 3. Here dP = A{p)6 {p^pi + m^) dP4^ is the volume element on 
the mass shell p'^pi = —m? in the momentum space. A{p) = 2, if is future directed 
and A{p) = otherwise; dP4 = y/^dp^dp^dp'^dp^. The quantity C[f] is Boltzmann's 
collision term. Its specific structure discussed e.g. by Ehlers will not be relevant for 
our considerations. Following Israel and Stewart we shall only require that (i) C is 
a local function of the distribution function, i.e., independent of derivatives of /, (ii) C 
is consistent with conservation of four-momentum and number of particles, and (iii) C 
yields a nonnegative expression for the entropy production and does not vanish unless / 
has the form of a local equilibrium distribution (see eq.(|^) below). Equation (|ID implies 
that the mass-m particles inbetween the collisions move according to the equations of 
motion 

where 7 is a parameter along the particle worldlines which for m > may be related 
to the proper time r by 7 = r/m. Since the particle four-momenta are normalized 
according to p^pi = — m^, the force has to satisfy the relation PiF* = 0. 

Both the collision integral C and the force F* describe interactions within the many- 
particle system. While C conventionally accounts for elastic binary collisions, we intend 
F* to model different kinds of interactions in a simple manner. Strictly speaking, F' 
should be calculated from the microscopic particle dynamics and, consequently, depend 
on the entire set of particle coordinates and momenta characterizing the A^-particle 
system. Since such kind of approach appears too complicated, in fact it already requires 
a solution of the many-particle dynamics, we follow here a different strategy. We assume 
F* to be an effective one-particle quantity which instead of depending on the coordinates 
and momenta of the remaining particles, is supposed to depend on macroscopic fluid 
quantities characterizing the A^-particle system as a whole. At the moment we do 
not specify this force. It will be determined below by general equilibrium conditions. 
Although lacking a detailed microscopic justification, the selfinteracting force concept 
will result in a transparent picture of the dynamics of the gas both macroscopically and 
microscopically. 



The particle number flow 4-vector N'^ and the energy momentum tensor T^^ are 
defined in a standard way (see, e.g., 0) as 

N' = I dPp'f {x,p) , T'' = J dPp^f {x,p) . (3) 

The integrals in (|^) and throughout the paper are integrals over the entire mass shell 
P^Pi = —m?- The entropy flow vector S'^ is given by 0, 

S^ = -Jp'^[f\nf-f]dP, (4) 

where we have restricted ourselves to the case of classical Maxwell-Boltzmann particles. 
Using well-known general relations (see, e.g., P]) we find 

N-a = j(c [f] - ^^'^) dP , T^t = Jp^(c [f] - mF^^)i dP , (5) 

and 

S-a = -j^^f{c[f]-mF^^)jdP. (6) 

In coUisional equilibrium, which we shall assume from now on. In / in (|^) is a linear 
combination of the collision invariants 1 and p". The corresponding equilibrium 
distribution function becomes (see, e.g., 0) 

/°(a;,p)=exp[a + /3„p'^] , (7) 

where a = a (x) and [3a (x) is timelike. Inserting the equilibrium function (|^) into eq.(|l|) 
one obtains 

P>a+Aa;;')PV = -"^Ai^' . (8) 

For a vanishing force P* the latter condition reduces to the "global " equilibrium 
condition of standard relativistic kinetic theory. For P* 7^ condition (H) is a 
"generalized " equilibrium condition (see fl^ and below). 
Use of Eq. (^ in the balances (|) yields 

Ar« = j F'fdP , T';^ = -m(3i J p^F'fdP . (9) 

For the entropy production density we find 

S-a = "^A / [a + In /°dP = -aiV.", - (3aT-.j, . (10) 

With / replaced by in (^) and (^), A^", and S"" may be split with respect to 
the unique 4-velocity u"" according to 

= nu" , T'^* = pu'^u'' + ph"'' , S" = nsu" , (11) 
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where h"-^ is the spatial projection tensor h"''' = g°'^ + m^u'', n is the particle number 
density, p is the energy density, p is the equilibrium pressure and s is the entropy per 
particle. The exact integral expressions for n, p, p and s are given by the formulae (177) 
- (180) in §. 

Using ( [TT| ) and defining 

r = / F7°dP , (12) 

n J 

the first eq.(|^) becomes 

h + Qn = riT , (13) 

where B = is the fluid expansion. It is obvious that F is the particle production 
rate. Similarly, with the decomposition (|11]) and the abbreviation 

t" = ml3i J p^F'fdP , (14) 

we obtain 



T^l + r = , (15) 



implying 



p+Q{p + p) =Uat\ {p + p)Ua + VaP = -hart' . (16) 

According to ([T3| ) and (|16D neither the particle number nor the energy-momentum 
are preserved in the presence of a nonvanishing F*. A varying particle number may 
be thought of as a result of inelastic interactions inside the medium |1^. Another 



possibility is particle production out of the gravitational field |jT9|, ^ |2T| . Source terms 
in the balances (|1^) are equivalent to the fact that the energy-momentum tensor T"-^ 
in (|ll]) is not conserved. Consequently, T"-^ is not the energy-momentum tensor for the 
selfinteracting g whole (see below). 

3. Quasilinear forces and generalized equilibrium 

The four-momenta may be decomposed into 

= Em" + Ae" , (17) 

where is a unit spatial vector, i.e., e'^e^ = 1, e'^Ua = 0. Consequently, one has 
E = —UaP'^ and A = CaP" and the mass shell condition p"'Pa = —m? is equivalent to 
= E"^ — m?. Moreover, hahP°'p^ = A^ is valid. For the force F^ we write analogously 
pm _ p ^x,p) -|- K {x,p) e"*, where F = —UmF^ and K = CmF^. The requirement 
PmF'^ = will be automatically fulfilled for 



E 
A 



F{x,p) . {U 
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With the famihar identification (3^ = Um/T it is obvious that only the part UmF"^ = —F 
contributes in the "source terms" ([T^) and ([T^). Up to now the function F{x,p) is 
completely general. Specific assumptions on the functional structure of F are necessary 
to evaluate explicitly the integrals in ( ]T^ ) and We will assume here that the 

macroscopic fluid dynamics in the presence of a nonvanishing force continues to 
imply only moments of the distribution function not higher than the second ones. It is 
obvious from the expressions ([T2|) and (0), which act as source terms in the balances 
(P!B|) and (|r^), respectively, that this requirement restricts F to depend on p°- only 
linearly. The most general linear expression for F is 

F {x, p) = Fo (x) + Fe {x) E + Fx{x)X. (19) 

Because of the additional dependence on through E = —UmP^ and A = CmP™ in the 
bracket in front of F (x, p) in (|18D, we call the force ([18| ) with the specific dependence (0) 
"quasilinear" . With this structure of the force the equilibrium condition (|^) becomes 

777 

P''a,a + Aa;;>)PV = Y [^0 {x) + Fe {x) {-Uap'') + F^ {x) (e.p")] . (20) 



Introducing the decomposition (|T^) also on the left-hand side of the last equation we 
find for the parameters a and (3 of the equilibrium distribution function 

Tfl Tfl 

a = -Fe , eVaa = -Fx , (21) 
where VaO = h^a^^^b, and 

P{a;b) = <P{x)gah , (22) 

with (j) = — Fo/(mT). The parameter a changes both in space and time while it is 
constant in the standard "global" equilibrium corresponding to the force-free case. With 
/3a = Ua/T the conformal Killing vector (CKV) property (p2D implies the relation 

0=^|, (23) 
which fixes the function Fq: 

777 

Fo = --e . (24) 

Fo vanishes for radiation (m = 0). 

It is remarkable that through Fq the force depends on the fluid expansion, a quantity 
characterizing the g whole on the macroscopic level. Similar dependences 

on macroscopic fluid quantities will be found below for Fe and Fx- Since both the 
microscopic particle four-momentum p°' and macroscopic fluid quantities, characterizing 
the system as a whole, enter the four-force, the latter represents a selfinteraction of the 
gas. 



The "sources" (p!2| ) and (|T^ in the balances (0) and (0), respectively, may be 
expressed in terms of Fq, Fe and F\. For V defined by (|12]) we find 



T^Z{Ef„ + f,^ (25) 



T \ n 

where M is the zeroth moment of the distribution function, M = J dP/° {x, p) = 
(p — 3p) /iTL^ . Analogously, in (14) may be written as 

777 

= -- [FoiV^ - FEUnT''^ + Fac^T'^"] , (26) 



implying 



TTl Tfl 

Uat" = - {nFo + pFe) , /leaf = -jpFxec . (27) 



The Gibbs equation 



Tds = d^+ pd- , (28) 
n n 

provides us with 

fi 

nTs = p-{p + p)- = Uat'' -{p + p)V (29) 
n 

for the time change s of the entropy per particle. Here we have used the balances (0) 
and dT^ ) to obtain the second equation (^) . 

With the identification a = p/T where p is the chemical potential and ||^ 

s='-^-^ (30) 
riT T ^ ^ 

the entropy production (|10D becomes 

S^^ = nVs + ns . (31) 

The equations of state for a gas may be generally written as (see, e.g., 0) 

p = nT , p = p{n,T) . (32) 

Differentiating the latter relation and using the balances (|T3D and (0) with (p9|) we 
obtain 

where the abbreviations 

dp _ (dp/dT)^ dp _ I dp\ 

'dp = (dp/dT)^ ' df=[df)^ 

have been used. 
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Restricting ourselves to the case s = 0, relation (|31| ) implies that there is entropy 



production only due to the enlargement of the phase space (cf 0, ^, 0). States of 
the system in which the gas particles (i) are governed by the equilibrium distribution 
function /°, (ii) the parameters a and Pa in /° obey (pOD, equivalent to (|21| ) and (|22|), 
and (iii) the entropy production 5." is given by 5." = nsT with F > 0, are called 
"generalized equilibrium" states |T^. We recall that the restriction s = refers to 



processes involving the force F*. Conventional dissipative processes due to deviations 
from collisional equilibrium have already been excluded by requiring the distribution 
function to have the structure (0) which makes the contribution of C [/] in the entropy 
production @ vanish. 

The condition s = relates the originally independent terms Uat"" and F: 

s = ^ Uaf" = {p + p)r . (34) 

Moreover, the CKV property (|2^) of Pi = Ui/T implies 

t e 

f = - 3 ■ P^' 
The last relation is consistent with (E^) for s = if 



r=fi-it^)e, (36) 



3 dp 



dp fmY ,_,,P±P_(P±P) . (37) 



where 



dp \T J nT \ nT 

Since ^ < 3 we have F > in an expanding universe, which is equivalent to 5"° > by 
virtue of (^TJ). The selfinteracting force provides us with a nonnegative expression for 
the entropy production in a similar sense in which Boltzmann's H theorem guarantees 
S°'a — /^'^ elastic binary collisions. Since relation (RBf) is a consequence of the 



equilibrium conditions it implies also the statement that except for p = p/3 (radiation) 
the generalized equilibrium conditions require a nonvanishing particle production rate F. 
If there is no particle production or if the rate at which the particle number changes is 
different from (|36|) , an equilibrium such as discussed here is impossible. 
The Gibbs-Duhem equation (see, e.g., 

dp=(p + p)^ + nTd(^) (38) 
together with the first equation of state ([3^ ) and equations ([T3|) and ( pS]) yields 




r - 1 - 7^ • (39) 
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Inserting here the particle production rate (0), the first equihbrium condition ( pTD fixes 
Fe: 



6m \p op J m \p op J m 



(40) 



For radiation, i.e. p = 3nT, the quantity Fe vanishes. For nonrelativistic matter with 
p = nm + ^nT, m ^ T one obtains 



(m > T) 



(41) 



As to the A part of the force Fx, the second equihbrium condition (pT]) aUows us to write 



T 



From the Gibbs-Duhem relation (|38D one finds that the spatial gradient Va (p/T) may 
be expressed by the fractional gradients of n and T, 

P\ pVaT 
p 

0, i.e., also the A-part of 



(43) 



,T/ n p T 

Since for radiation n oc is valid, one has Va{p/T) 
the force vanishes in this limit. The motion of ultrarelativistic particles is force free. 
This corresponds to the vanishing of the source terms (|25| ) and (|26| ) in the macroscopic 
balances ( [1^ ) and (|15|), respectively. 
For nonrelativistic matter we obtain 

"KTJ T T 
and the A part of the force becomes 



(m >T) 



(44) 



2^2 



(m > T) 



(45) 



where the relation VaT/T + iia = 0, which follows from the CKV property (p^), has 
been used. 

Combining (g^), (|0D and (g^), we find for the quantity F = -UmF"^ in (|1|) 

dp\e 



F{x,p) 



3 m 



P 



p 



p 



p"" . (46) 



Its first two parts are proportional to the fiuid expansion while the last part is 
proportional to the spatial gradient of p/T. While for ultrarelativistic particles 
(p = 3nT) each of the three parts in (^6]) vanishes separately, the force is different 
from zero for any other equation of state. We conclude that a system of particles with 
nonzero mass m is in generalized equilibrium if inbetween the elastic, binary collisions 
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the particles move in the force field ( fZ^ with (^4^)- the nonrelativistic limit m ^ T 
the expression (^61) reduces to 



F {x, p) = {E-m)^ + Wu^'iia {m:^T) . (47) 

We recall that this force is not an external force but determined by the macroscopic fluid 
quantities of the system of colliding gas particles themselves. The force F"^ represents 
a selfinteraction of the gas under the equilibrium conditions (0) and (|22|) . 

The discussion so far allows us to reconsider earlier work ( 0, H, |1^ ) on kinetic theory 
and particle production based on a Boltzmann equation modified by a "source" term 
which was assumed to be linear in the distribution function with an additional linear 
dependence on the particle four-momenta of the factor in front of /° (see formulae (37) 
and (38) in 0]). While the second assumption is equivalent to our eq. (|T9|), there 
is no need for a conterpart of the first one in the present paper. Although formally 
similar, the description in terms of selfinteracting forces is physically more transparent 
than the somewhat vague concept of a "source term" in Boltzmann's equation. For 
comparison, we list up the correspondences between the spacetime functions z/, r and 
a, characterizing the mentioned source term in ||, [l^, and the quantities Fq, Fe and 
Fx of the present paper: 

m„ 1 m ^ 1 m ^ , ^ . 

ly^-Fo, -^TfFE, -^TfFx. 48 

The "sources" z/, 1/r and 1/cr (see 0, ||, |T^) may be understood in terms of 
selfinteracting forces on the microscopic gas particles. 

The main advantage of the force concept compared with the "source term" approach 
is the following. Once determined by the equilibrium properties of the gas, the force 
may now be used to reconsider the microscopic equations of motion. With the splitting 
(p!7|) the left-hand side of the second equation (^ may be written as 

Dp' dE , Du' dX , Be' 
-^M + E— h -re + A- 



dr dr dr dr dr 

Contraction with Ui yields 

Dp' dE ^ De' dE ^ -Dm 
Ui— — = — - — h Xui— — = — Ae 



Use of 



dr dr dr dr dr 



Du' . 
dr ' m 



allows us to write 
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Applying here the decomposition |2^ 



Ui.,n - -UiUn + + UJ^n + , (49) 

where ujab = h^h'^uic-iti and aat = hlhfu{^c;d) — jhab, the equation of motion (H), projected 
in direction of m*, 

= ^^F' = -F, (50) 
may generally be written as 

^ + —e^u. + -eVV.„ + fe = F. (51) 
(IT m m Sm 

In the following section we will integrate this equation for specific cases. 



4. Fluid and particle dynamics 

Obviously, all relations of the preceeding section are applicable on an arbitrary, 
given gravitational background. Considering the gravitational field of such kind of 
selfinteracting equilibrium system requires further investigations. Since the energy 
momentum tensor T"* in (|^) and ([TI|) is not conserved, it is not a suitable quantity 
on the right-hand side of Einstein's field equations. The question arises whether 
one may introduce an effectively conserved energy-momentum tensor T""^ satisfying 
= T"^ + t°' = 0, where the "source" is mapped onto suitable components of the 
energy-momentum tensor T°'^. This probem has been solved previously P, ^ [1^ and 
the corresponding results may be applied immediately. The essential point is that the 
balances ( |T^ ) following from (p!5D are identical to the balances 

p + e(p + p + 7r) = , {p + p + Tr)Ua + Va{p + Tr) = (52) 
following from T"^ = with 

f"'' = pu^u^ + (p + tt) h"^ (53) 
and the identifications 

Uat" = -OtT , haif = Tlila + VaTT . (54) 

The source terms in ([15|) may consistently be mapped onto an effective viscous 
pressure vr of a locally conserved energy- momentum tensor T°*. This viscous pressure 
is determined by eqs. (B3) and (E^): 
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The quantity T"'' in (^3|) with ( |55D may be regarded as the effective energy-momentum 
tensor of the selfinteracting many-particle system. 

The first purpose of this section is to demonstrate how vr is determined by the 
selfinteracting force. Combining ( PTj ) with (|3^) and ( |5BD yields 

{nFo + pFE)e-' . 



TT = — 



Applying (|3^) with 
for TT are possible: 



(56) 



and (|40|) as well as (0), the following equivalent representations 



IT 



m . , 
^ (P + P) 



M ■ 
— Fo + Fe 
L n 



e 



-1 



my 

TJ 



p-3p p 



1 



(57) 



nm nm 

While (|36D and (|55|) are previous results [H, ||, the expressions (|^), (|56|) and (^) 
are new. The latter are equivalent to (|3^) and (^) but differently from these previously 
obtained formulae the relations (p5|) and (|56|) (or (|57|)) establish a direct connection 
between the force on the gas particles and the particle production rate F and the 
effective viscous pressure vr, respectively. The deviations from the standard perfect fluid 
behaviour, represented by vr, are immediately connected to a microscopic selfinteraction 
of the gas, which governs the particle motion according to (|5lD. 

The backreaction of a nonvanishing vr on the fluid dynamics ([5^ ) depends on the 
equations of state. There is no backreaction for p = 3nT. But for any other equation of 
state the force F"* and, consequently, the particle production rate T and the effective 
viscous pressure vr, are different from zero. The backreaction is largest for m ^ T, 
i.e., for nonrelativistic matter. Introducing a length scale a according to = 3d/a one 
obtains 

r=®, ^ = -f' ' P ^ (m > T) . (58) 

The temperature changes according to (pSf ) for any equation of state. The familiar 
behaviour of nonrelativistic matter characterized by n oc a~^, T oc and p oc 
for r = TT = is modified considerably. All thermodynamic quantities decrease more 
slowly for TT 7^ since the decay of n, T and p due to the expansion is counteracted 
by corresponding production terms reflecting the selfinteraction of the gas particles on 
a macroscopic level. Assuming that matter in this kind of equilibrium dominates the 
cosmological dynamics, the expansion of the universe is modified as well [§, In a 



homogeneous and isotropic universe the lenght scale a coincides with the scale factor of 
the Robertson- Walker metric and, in the spatially flat case, obeys the equation 

a2 



3— = , 

where k is Einstein's gravitational constant. For radiation with vr 
course, a oc t^/^. Inserting, however, the energy density p from Eq. 



(59) 

we recover, of 
into Eq. (g) 
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we find tliat tlie scale factor behaves such as 

aoct^/^ (m>r) , (60) 

i.e., a > instead of the familiar a oc t"^^^ with a < for p oc a^^ corresponding to 
r = 0. Generalized equilibrium of a selfinteracting massive gas universe implies power 
law inflation. The possibility of an accelerated expansion of the universe may be traced 
back to specific forces on the microscopic constituents of the cosmic medium. 

Since these forces are explicitly known, we may solve the microscopic particle 
dynamics as well. For homogeneous and isotropic universes, characterized by iii = 
ain = 0, equation (|5ll) reduces to 

With dr = dt (m/E), \^ = E"^ — m^, and dE/dt = E, the last equation is equivalent to 

2~ ~ ~ 2^ ■ ^^'^l 

E'^ — a -C/ — 

In the limit F = (geodesic motion) we obtain 

E^-m' = \^ (X , (F = 0) , (63) 

implying the expected behaviour E oc for massless particles (photons) while the 
nonrelativistic energy e = E — m with e <^ m of massive particles decays as e oc a~^. 
The most interesting case ( ^7D simplifies to F = [E — m) Q/3 for iii = 0. Under these 
conditions which macroscopically are characterized by (|58D and (|60|), the solution of 
equation (pT| ) is 

F - m oc a"\ (^F = (F - m) ® ^ . (64) 

The nonrelativistic energy of massive particles under generalized equilibrium conditions 
in a (quasi-) linear force field decays linearly with the cosmic scale factor, i.e., the 
selfinteracting force makes nonrelativistic particles behave like radiation. Massive 
particles under the action of a force field ^T^ ) with F = (F — m) 6/3 behave like photons. 
This is the microscopic counterpart of the statement that radiation and nonrelativistic 
matter may be in equilibrium in the expanding universe, provided the number of 
matter particles increases at a specific rate P]. The solution (|6^) for the microscopic 
particle dynamics together with the corresponding macroscopic fluid behaviour (^) 
and ( |60D constitutes a nontrivial, exactly solvable model of a gas universe in generalized 
equilibrium. 



5. Summary 



We have investigated the kinetic theory of a classical relativistic gas of particles moving 
under the influence of a simple, selfinteracting force inbetween equilibrium establishing 
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elastic, binary collisions. The concept of a selfinteracting force was shown to provide a 
comprehensive picture of the gas dynamics both macroscopically and microscopically. 
On the macroscopic level it gives rise to an effective viscous pressure which, on the 
microscopic level, corresponds to a deviation from geodesic particle motion. A deeper 
understanding of previously obtained results on the backreaction of particle production 
processes on the cosmological dynamics was obtained. We established an exactly solv- 
able model according to which massive gas particles under the action of a specific force 
behave hke photons. A universe of gas particles of this kind is necessarily in a state of 
accelerated expansion corresponding to a behaviour a oc t^/^ of the scale factor a of the 
Robertson- Walker metric. 
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